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We develop a theory for spin transported by coherent Neel dynamics through an antiferromagnetic insulator 
coupled to a ferromagnetic insulator on one side and a current-carrying normal metal with strong spin-orbit cou¬ 
pling on the other. The ferromagnet is considered within the mono-domain limit and we assume its coupling to 
the local antiferromagnet Neel order at the ferromagnet|antiferromagnet interface through exchange coupling. 
Coupling between the charge current and the local Neel order at the other interface is described using spin 
Hall phenomenology. Spin transport through the antiferromagnet, assumed to possess an easy-axis magnetic 
anisotropy, is solved within the adiabatic approximation and the effect of spin current flowing into the ferromag¬ 
net on its resonance linewidth is evaluated. Onsager reciprocity is used to evaluate the inverse spin Hall voltage 
generated across the metal by a dynamic ferromagnet as a function the antiferromagnet thickness. 

PACS numbers: 72.25.Mk, 75.47.-m, 75.50.Ee, 76.50.+g 


Spintronics of antiferromagnets (AFs), where AFs take on 
the role of the central active component, is identified as one of 
the most important emerging topics in the held of magnetism 
today m. Robustness to magnetic perturbations due to their 
total magnetic compensation, as well as characteristic dynam¬ 
ical scale in the THz range may render AFs advantageous over 
ferromagnets (Fs) for spintronics device applications. In ad¬ 
dition, recent works on AFs have shown that the important 
phenomena responsible for the success of F-based spintron¬ 
ics also have AF counterparts, giving added impetus for AF- 
based spintronics research. Indeed, giant magnetoresistance 
and current-induced torques j2[, anisotropic magnetoresis¬ 
tance 0 and spin superfluidity J4], as well as current-induced 
domain wall motion Q and coupled dynamics between con¬ 
duction electrons and background magnetic texture 0, are all 
shown to be possible in AFs as well. 

An important aspect of AF-based spintronics is the use of 
AFs as a medium to transport spin angular momentum. Spin 
transfer through AFs has been the focus of several recent ex¬ 
perimental endeavors. Both Hahn et al. j7| and Wang et 
al. 0 demonstrated spin transport through an AF insulator, 
NiO, using an YIG|NiO|Pt heterostructure (YIG standing for 
the insulating ferrimagnet yttrium iron garnet). Inverse spin 
Hall signal showed robust spin pumping from YIG into Pt 
even in the presence of the intervening NiO, suggesting effi¬ 
cient spin transport through the AF. More recently, Moriyama 
et al. used spin-torque ferromagnetic resonance (ST-FMR) 
to demonstrate the propagation of spin excitations through a 
metallic AF, IrMn, using a Pt|IrMn|CoFeB trilayer 0 as well 
as NiO using a Pt|NiO|FeNi trilayer Col. Spin current injected 
from the Pt was shown to change the FMR linewidth, also sug¬ 
gesting the transfer of spin angular momentum through the 
central AF. Given the rising interest in AF spintronics and 
the recent experimental focus, a theoretical account of spin 
transport through an experimentally relevant normal metal 
(N)|AF|F trilayer is highly desirable. 

In this Letter, we develop a general phenomenology for spin 
transport through an AF by collective Neel order parameter 



effective junction 


FIG. 1. Normal-metal (N)|antiferromagnet (AF)|ferromagnet (F) tri- 
layer considered in this work. N sustains a dc charge current j and 
F is described by a time-dependent macrospin S(t). Spin transfer 
( J‘ x ) occurs via the exchange coupling J at the AF|F interface, while 
spin transfer across the AF|N interface has a spin transfer torque con¬ 
tribution ( J* tt ) (proportional to the effective interfacial spin Hall an¬ 
gle i?) and a spin pumping contribution (J* p ) (proportional to the 
interfacial spin-mixing conductance cr-ly The AF Gilbert damping, 
parametrized by a', leads to the loss of spin current (J^.) in the AF 
bulk. The central AF can be thought of as an effective interface that 
couples j and S with an effective spin Hall angle iT. 


dynamics, focusing on an N|AF|F trilayer relevant for both the 
spin-pumping/inverse spin Hall as well as the ST-FMR exper¬ 
iments mentioned above (see Fig. |T]». Spin Hall phenomenol¬ 
ogy, applicable to a wide range of different AF|N interfaces 
obeying certain structural/crystalline symmetries, is utilized 
to model the spin transfer at the AF|N interface, while the ex¬ 
change coupling is assumed at the AF|F interface. As one of 
the main achievements of this work we develop a simple “cir¬ 
cuit” model, a pictorial visualization of spin flow, that allows 
one to keep track of spin transfer through various parts of the 
heterostructure (see bottom half of Fig. [TJ. From the circuit 
model, we see that spin is both injected into (i.e., (J s ' |t )) and 















2 


ejected (i.e., (JL)) out of the AF at the AF|N interface due to 
spin Hall/spin-torque effects and spin pumping, respectively. 
The collective Neel dynamics leads to Gilbert damping and to 
the loss of spin current (i.e., within the AF bulk, and the 
exchange coupling at the AF|F interface leads to spin transfer 
(i.e., (J e 5 x )) across the interface. We first study how spin trans¬ 
port through the AF modifies the linewidth at FMR, akin to the 
ST-FMR MM- The FMR linewidth is quantified in terms of 
the effective spin Hall angle and spin-mixing conductance at 
the AF|N interface, the exchange coupling at the AF|F inter¬ 
face, as well as AF Gilbert damping. We show that linewidths, 
measured for different electrical currents in N and AF thick¬ 
nesses, can be used to extract the effective spin Hall angle 
and spin-mixing conductance at the N|AF interface, as well as 
the bulk Gilbert damping. By invoking Onsager reciprocity, 
we also make connections with the inverse spin Hall experi¬ 
ments 00 and compute the inverse spin Hall voltage gener¬ 
ated across N arising as a result of a dynamic F macrospin. 

As shown in Fig. [T] an insulating AF is attached on one side 
to a mono-domain F and on the other to a paramagnetic N with 
strong spin-orbit coupling. The N and F sustain dc charge 
current density j and a time-dependent macrospin S(t), re¬ 
spectively. We consider a bipartite AF, which can be charac¬ 
terized by two hydrodynamic variables, n(x, t) and mix, f), 
parametrizing the staggered (Neel) and smooth (magnetic) 
components of the spins, respectively 03]. We assume easy- 
axis magnetic anisotropy along the z axis in the AF, as well 
as full translational and rotational symmetries in the yz planes 
so that our treatment essentially reduces to a one-dimensional 
problem that depends only on the coordinate x. The free en¬ 
ergy i? for the AF and its coupling to the F reads 


& = 



[d x n(x)l 2 
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K 7 ) 
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- JS ■ n(L ), 
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where A and x are the Neel order stiffness and spin suscepti¬ 
bility, respectively, k the uniaxial anisotropy parameter, and J 
the exchange coupling between AF and F. 

The Landau-Lifshitz-Gilbert dynamics in the bulk AF cor¬ 
responding to Eq. 0 can be written as 

s(n + anx rh) = x'mx n, ( 2 ) 

s{rn + am x rh + a'n xn) — nx (Ad 2 x n + Kn z e z ), (3) 

where a and a' are (independent) Gilbert damping param¬ 
eters and s is the roughly saturated spin density (per unit 
length) 02- In the limit of slow dynamics (i.e., r » ti/E ex , 
where r is the time scale for AF dynamics and E ex is the 
AF exchange energy scale) and strong local Neel order (i.e., 
\mix, f)| <s 1), one may first use Eq. 0 to solve for m, insert 
the solution into Eq. 0, and arrive at the following dynamics 
for the AF Neel vector, 

Xs 2 n x n + sa'n x n « An x d 2 n + kux n z e z . (4) 


The AF spins are excited by the current and the dynamic 
macrospin. The effects of these external perturbations are lo¬ 
calized at the interfaces and thus enter the AF dynamics as 
boundary conditions. At the AF|N interface, spin currents 
arise via spin transfer torque and spin-pumping, both of which 
can be accounted for using spin Hall phenomenology lfT3l . 
Based on structural symmetries at the interface, spin Hall phe¬ 
nomenology allows us to write down a general expression for 
the spin transfer torque that applies to a variety of F- and AF- 
based heterostructures with different microscopic details. In 
the presence of full translational and rotational symmetries 
in the yz plane and with the breaking of reflection symmetry 
along the x axis, there are two contributions to the spin current 
(integrated over the interface area) flowing into AF m3 

Ji = [&n x (e x x j)xn- ha n n x n]\ x=0 = J s „ - J sp , (5) 

where the first term is the so-called spin Hall-like (dissipative) 
contribution and the second term describes spin-pumping. 
The coefficient 8 is proportional to the (tangent of the) ef¬ 
fective spin Hall angle at the AF|N interface 03; although i? 
can, in general, depend on the orientation of n, we will treat 
it as a constant here. We will also disregard any anisotropies 
of ar^i w ith respect to the orientations of n and n, assuming 
that the exchange energy scale at the interface dominates over 
the energy scale of spin-orbit interactions. 

While the coefficient a ^ is proportional to the real part of 
the (generally complex) spin-mixing conductance g 1 ' J - for the 
AF|N interface |T3 . its imaginary part gives rise to a term 
oc ft in Eq. 0. The assumed structural symmetry, in prin¬ 
ciple, also allows for the so-called field-like (reactive) con¬ 
tribution given by rj(e x x j ) x n. Assuming the magnitudes 
of all four terms are small (as compared to yAx to be pre¬ 
cise) we will compute the FMR linewidth to linear-order in 
the coefficients f), i) and Since the terms that do not ap¬ 
pear in Eq. ([5]) are non-dissipative terms, they do not enter the 
linewidth at linear-order in these coefficients. We will, there¬ 
fore, drop these contributions from the following discussion 
and only consider the two terms in Eq. ([5]). 

The exchange coupling J [given by the last term in Eq.(|T]»] 
leads to the following expression for spin current flowing 
across the AF|F interface, 

8 S? 

J r^ J L=-ggXS = J Sxn\ x=L . ( 6 ) 

The spin current in the AF bulk can be read off from Eq. Q 
(dropping the Gilbert damping and anisotropy terms) and the 
resulting continuity equation [i.e., srh = -d x {-Anxd x n)], so 
that J^ F (x) = -A n(x) xd x n(x). The continuity of spin current 
across each interface then leads to the boundary conditions 

Ji=J S A f(* = 0), J; = J1f(x = L). (7) 

The dynamic Neel texture n(x, t ) can be obtained using the 
low-frequency (adiabatic) approximation, valid in the regime 
Q <s Qy, where £1 and Go are the FMR and the AF resonance 
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FIG. 2. The interfacial contribution Sa^ (dashed lines), the bulk 
contribution <5or® (dotted lines) and the total contribution Sa F (solid 
lines) to the extrinsic FMR linewidth (S set to unity) are plotted as a 
function of the (normalized) system size L/A. We fix the following 
parameters: r] = 1, ftj y /bo = 0.01 and a = 0.01. Two regimes 
are considered for the AF Gilbert damping a: (a) the strong damping 
regime a = 0.2; and (b) the weak damping regime a = 0.01 (see text 
for more details). 


frequencies, respectively. Within this approximation, the AF 
Neel texture is first solved for an arbitrary static S, the result 
denoted by n (0) (x, S). Since Sit) varies sufficiently slowly in 
time compared to the characteristic AF time scale, the Neel 
texture in the adiabatic limit will arrange itself into the static 
configuration corresponding to Sit ) at every moment in time 
and is given by n(x, t ) as n (0) [x, S(t)] = n (0> (x, t). The above 
calculation does not account for spin current losses due to the 
AF dynamics (i.e., spin-pumping at the AF|N interface and 
Gilbert damping in the AF bulk). Taking these losses into ac¬ 
count up to linear-order corrections to the adiabatic result, the 
spin current (J e ' x ) entering F, time-averaged over a cycle of 
FMR precession (the angle brackets (■ • ■) hereafter represent¬ 
ing time-average over a cycle of FMR precession), is given 

b y ( J L) = m - <J S p) - < j g> (c - f - Fi g-0’ where the s P in ~ 

transfer torque contribution is given by inserting the adiabatic 
result for the Neel texture into Eq. <0 

<J4> = 7<S(f)xn (0) (L,f)>, (8) 

and the loss terms read 


tained if we consider small deviations of S(t) away from 
the z axis (parallel to the static FMR field and the AF easy- 
axis); we take j = j y e y + j z e z and assume \j\ to be weak 
such that a linear-response treatment is sufficient. In this 
case, the Neel unit vector n should not deviate far from the 
z axis and we may evaluate the above results with respect 
to small transverse fluctuations, i.e., S(t ) ~ S [s x (f), s y (t), 1] 
and n(x,t ) & [n r (x,t),n y (x,t),l] with |.S’ l (f)|, |.v y (t)| <K 1 and 
| n x (x,t)\,\n y (x,t)\ «: 1. Within this treatment, the trans¬ 
verse components n ± = (n x , n y ) T obey ^s 2 fi ± + san ± = 
Ad 2 x n ± - Kn ± [c.f. Eq. (j4|], and n i0) (x, t) has the form 

n (0) « e, + f(x)[e z xS(t)]xe. + g(x)e z x S(t) + h(x)e x , (10) 

where the functions f(x) and g(x) (to linear-order in the cur¬ 
rent) are given by 


f(x) = 


1 cosh | 

S cosh - sinh ’ 

A 7] A 


1 sinh ^ + jj cosh ^ {fj y 
S (cosh \ + i sinh \) 2 


( 11 ) 

( 12 ) 


and li(x ) oc § j z is not explicitly shown here since this term will 
not contribute to the linewidth within the current theoretical 
treatment. Here, jj = JS / VAk, and A = \JA / k is the AF 
healing length. 

The spin current J e ' x entering F modifies the F dynamics as 


hS = bx S - T ^-S x S + J°, 
S 


(13) 


where a F is the intrinsic Gilbert damping parameter in F and 
b = -b( t e : is the static FMR field (in units of energy). Insert¬ 
ing Eq. m i nt0 Eqs. (|8j and (|9]» and performing the time- 
average over the last two terms in Eq. ( fT3] >, the full FMR 
linewidth can be read off directly by summing the coefficients 
appearing in front of ( S x S ). The total Gilbert damping pa¬ 
rameter is then given by a' F = a F + da ( F + 6a^ } = a F + 5a F , 
where the extrinsic contribution day has the interfacial contri¬ 
bution 6a^ and the AF bulk contribution da 


da { F = 
6a { F = 


■ 

^ (cosh j + - sinh j) 
a \ + \ sinh f 
$ (cosh j + i sinh 


(14) 

(15) 


< J/ p ) = ha n (n <0 \ 0, t) x n (0) (0,0), 

r L (9) 

( Jq) = sa' I dx (n (0) (x, t) x fi {0 \x, t)). 

Jo 

The first term in Eq. (|9]> describes (time-averaged) spin current 
lost due to spin pumping at the AF|N interface and the second 
term corresponds to Gilbert damping in the AF bulk. 

An analytical result for the FMR linewidth can be ob- 


where a = sa'A/2%. The former originates from spin injection 
and spin-pumping at the AF|N interface while the latter from 
Gilbert damping in the AF bulk. Eqs. ( |T4| and ( fl5] > constitute 
the main result of this work. 

As seen from Eqs. ( [T4| and ( p~5| >, the healing length A sets 
the distance over which spin propagation decays inside the 
AF. The healing length is determined from the slope of the 
linewidth vs. j y curves for various thicknesses L and by ex- 
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trading the decay length. It is important to note that the cur¬ 
rent theory only considers spin transport mediated by coherent 
Neel dynamics, and does not take account of spin transported 
by incoherent thermal magnons. The latter contribution is 
suppressed at sufficiently low temperatures by some power of 
the ratio T/T N , where TV is the Neel ordering temperature 
of the AF, and since magnon-mediated transport is expected 
to decay over the spin diffusion length /f sc j, it is strongly sup¬ 
pressed for d s( i <sc L. 

Once the AF healing length is known, the parameters ft, ar^l 
and or can be extracted by measuring the FMR linewidth for 
various j y and L. While the effective spin Hall angle ft can be 
obtained from the slope of a linewidth vs. j y curve, the Gilbert 
damping parameter a can be extracted in the regime L » A, 
in which the linewidth depends only on a (see Fig. [2|, 

(16) 


For L «: A, we expand 5a F to linear order in L/A, 


6a F ~ 7 T (+ Or U 


S\b 0 


= c 0 + Ci~ , 
A 


a - 

1 


2(ftj 
bo 


y 11 
■ +a li 


(17) 


from which we see that can be extracted at zero current 
(i.e., j y = 0) and measuring the linewidth for L <sc A. 

The extrinsic linewidth exhibits qualitatively different be¬ 
havior depending on the relative magnitudes of the bulk and 
the interfacial contributions (see Fig. |2j. For strong Gilbert 
damping (ci > 0) the bulk damping in the AF dominates over 
the interface effects and the linewidth grows initially as L in¬ 
creases, saturating eventually as L/A —> oo [see Fig. [2ja)]_ In 
the limit of weak Gilbert damping [see Fig.|2]b)], i.e., ci < 0 
(and 6ap < c o), 6a F exponentially decays as L increases. In 
Ref. [8) the FMR linewidth was measured for various AF thick¬ 
nesses in the absence of the electrical current. The gradual 
increase in the linewidth obtained there as a function of the 
thickness is more consistent with our strong Gilbert damping 
regime [c.f. Fig.[2]V)]. 

Moriyama et al. ITOl has recently reported an ST-FMR 
measurement for single crystal NiO using a Pt|NiO|FeNi tri¬ 
layer, observed a linear dependence of the FMR linewidth 
on the Pt electrical current, and extracted the slope [c.f. 
Eq. {H}]. The reported ratio fi ~ 0.82 between the slopes 
for the Pt|NiO|FeNi and Pt|FeNi samples indicated relatively 
efficient spin transfer through NiO compared to previous re¬ 
ports 01 , and was attributed to the single crystal nature of 
the AF layer. Applying Eq. ( [T4| directly to the experiment 
(see also Ref. EB, the theoretical ratio J3 t h (ignoring a pos¬ 
sible difference in ft between the two samples) is given by 
/T t h = [cosh(L//i) + sinh(L//l)/77] -2 . Using A ~ 100 nm, q ~ 1 
(appropriate for NiO and NiO|FeNi interface) and L ~ 10 nm, 
we obtain /J,h ~ 0.8, which is quantitatively consistent with 
the experiment fTOl . 


Our results can be used to make a connection with the recip¬ 
rocal experiments 00 , in which spin transfer through the AF 
is quantified by measuring the inverse spin Hall voltage Vishe 
generated across N by a dynamic F (see Fig. [TJ. From spin 
Hall phenomenology and Onsager reciprocity, the electromo¬ 
tive force generated in N is given by e = ft(nxn)xe x \ x =o El- 
Utilizing the adiabatic result n (0 \x,t) [i.e., Eq. ©] with 
j = 0, the (time-averaged) motive force becomes (e) = 
-(ft*6 2 b 0 /H)e y , where 6 « (s 2 + s 2 ) 1/2 is the cone angle, 
ft* = #/{yV[ C0S h(T//l) + sinh(L//l)/77] 2 }, and jj- is the vol¬ 
ume of F. This leads to an inverse spin Hall voltage 


&e 2 b 0 { 

Vishe =---- , 

ny F (cosh \ + i sinh 


(18) 


where l is the length of N in the y direction. We can arrive at 
the same result by treating the central AF as an effective junc¬ 
tion between the N and F subsystems (see Fig. [T]). Namely, 
from Eq. the macroscopic coupling between current j 
and the macrospin dynamics in F is given through the torque 
t = i TO 2 j y e z + (term oc j.) acting on the latter, where ft* is 
the overall torque coefficient for the effective junction. By On¬ 
sager reciprocity, this torque gives rise to the inverse spin Hall 
voltage Eq. ©■ Within the current theory, Eq. indicates 
that the decay length for Vishe as the AF thickness increases 
is set by A. 

We note in conclusion that the current work considers spin 
transfer purely mediated by the coherent Neel dynamics, cor¬ 
responding to the so-called superfluid contribution to spin 
transport. As the relevant experiments are performed at room 
temperature 0Ema reconsidering AF spin transport by ac¬ 
counting for the incoherent thermal magnons and studying 
their effect on the FMR linewidth would be valuable, and 
will contribute to the general understanding of the “two-fluid” 
(condensate and thermal cloud with mutual interactions be¬ 
tween them) nature of spin transport via collective excitations 
in an AF. 
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